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Continuous-variable quantum process tomography with squeezed-state probes 
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We propose a procedure for tomographic characterization of continuous variable quantum opera¬ 
tions which employs homodyne detection and single-mode squeezed probe states with a fixed degree 
of squeezing and anti-squeezing and a variable displacement and orientation of squeezing ellipse. 

Density matrix elements of a quantum process matrix in Fock basis can be estimated by averag¬ 
ing well behaved pattern functions over the homodyne data. We show that this approach can be 
straightforwardly extended to characterization of quantum measurement devices. The probe states 
can be mixed, which makes the proposed procedure feasible with current technology. 

PACS numbers: 03.65.Wj 


I. INTRODUCTION 


As the complexity of quantum information processing 
devices increases, there is a growing need for tools for 
their characterization and benchmarking. Quantum op¬ 
erations and channels can be completely characterized by 
quantum process tomography m, which represents an 
extension of quantum state tomography IHS] to quan¬ 
tum operations. Typically, the quantum operation £ is 
probed with a sufficient number of input states pj , mea¬ 
surements in several different bases are performed on the 
output states £ {pj ), and the quantum operation is recon¬ 
structed from the experimental data. Alternatively, in 
the ancilla-assisted quantum process tomography EHn] 
the operation £ is probed with one part of a single hxed 
entangled bipartite state pab and £ is determined from 
measurements on the output bipartite state. This latter 
approach is based on the Choi-Jamiolkowski isomorphism 
[HIS], which tells us that if the probe state is pure and 
maximally entangled, |$)ab = then 

the output bipartite state xab = 1a ® £’b(|$)($|) is 
directly isomorphic to the operation £. Here d denotes 
the dimension of input Hilbert space 'Hin and the states 
\j) form an orthonormal basis in "Hin- 

Quantum process tomography works particularly well 
for few-qubit systems, and it has been successfully ap¬ 
plied in the past to characterization of various single¬ 
qubit and two-qubit operations [TTl [T^ [TBH^ . As the 
number of qubits N increases, the full tomography be¬ 
comes challenging since the number of parameters that 
have to be estimated grows exponentially with N. In 
some cases, scalable quantum process reconstruction may 
be achieved e.g. by approximating the operator x by 
a matrix product state [12 [Hj or by using compressed 
sensing techniques [24U26] . 

Besides the issue of Hilbert space dimension, the quan¬ 
tum process tomography is also affected by the range of 
practically accessible input probe states. This is partic¬ 
ularly relevant for continuous variable quantum process 
tomography 0 0 HTHSH], which aims at characteriza¬ 
tion of quantum operations on modes of quantized elec¬ 
tromagnetic fields. Here, the most natural and read¬ 
ily available probe states are represented by coherent 


states |a), and the output states can be conveniently 
measured with homodyne detectors 00. Recently, this 
approach has been successfully employed to character¬ 
ize a single-mode lossy channel m, and a conditional 
single-photon addition and subtraction [30]. Moreover, 
the coherent states were also used as probes for com¬ 
plete tomographic characterization of single-photon de¬ 
tectors jSTHSS]. Probing quantum processes with coher¬ 
ent states essentially amounts to determining a Husimi 
Q-function of the operator x- More precisely, assuming 
that the measurements on output states are described by 
a POVM with elements Hj, the probability of measure¬ 
ment outcome H^- for input probe coherent state la) reads 
Pj{a) = Tr[|a*)(a*| 0 Hj y]. Usually, one would like to 
reconstruct the matrix elements of x in Pock basis. To see 
the connection between the Husimi Q-function and the 
matrix elements in Fock basis, recall that the Q-function 
of an operator A is defined as Q{a) = (a|A|a)/7r. We 
have 


-|a|^ ^ a*™a" 

Q{a,a*) = - - X 


( 1 ) 


which shows that the Q-function is a generating function 
of matrix elements Ajn,n = {m\A\n) in Fock basis [28], 


A — 


^m+n 




Q{a, a*)e 


|a|^ 


• ( 2 ) 


Q; = Q ;*=0 


Here a and a* are formally treated as independent vari¬ 
ables. Estimation of Am n from experimental data re¬ 
quires inversion of Eq. Q when the Q-function is not 
known precisely. This can be a delicate procedure sensi¬ 
tive to statistical fluctuations of the data. 

In Ref. [27], elements of quantum process matrix of 
a lossy channel were reconstructed from the experimen¬ 
tal data with the help of regularized version of Glauber- 
Sudarshan P-functions [M] of operators |TO)(n|. This 
approach approximates the calculation of derivatives in 
Eq. by evaluation of a suitable linear combination of 
the experimental data. In Ref. [31], POVM elements of 
a single-photon detector were reconstructed from mea¬ 
surements on probe coherent states by solving a convex 
optimization problem that included an extra constraint 
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which ensured a smooth structure of the reconstructed 
POVM elements. Later on, maximum-likelihood estima¬ 
tion was employed for reconstruction of quantum opera¬ 
tions and detectors probed with coherent states [ 301133 ] ■ 
This latter approach avoids the complications with di¬ 
rect linear inversion ([^, but it requires some truncation 
of the infinite-dimensional operator y. 

In this paper, we investigate characterization of con¬ 
tinuous variable quantum operations which is based on 
single-mode squeezed probe states and homodyne detec¬ 
tion on output states. By using squeezed states instead 
of coherent states we avoid the problems with linear in¬ 
version of the data and we show that the matrix elements 
of quantum process x in Fock basis can be determined 
by averaging suitable well behaved pattern functions [ 371 - 
I39| over the homodyne data. Our procedure assumes that 
all probe states have the same variances of squeezed and 
anti-squeezed quadratures, and these variances need to 
be known and kept constant during the whole measure¬ 
ment. The probe states also need to be phase shifted and 
coherently displaced in a controlled way, which is feasi¬ 
ble with current technology. Importantly, our procedure 
works for realistic mixed squeezed states and the only re¬ 
quirement is that the variance of the squeezed quadrature 
is below the coherent state level. Although we focus on 
linear reconstruction based on the formalism of pattern 
functions, the data could be processed by other means, 
such as the maximum likelihood estimation. Our work 
provides an important insight into the utility of squeezed 
states for tomography of quantum processes. 


II. QUANTUM PROCESS TOMOGRAPHY 

In what follows we shall consider chracterization of 
a single-mode quantum operation £. According to the 
Choi-Jamiolkowski isomorphism |14[ 115] . such operation 
can be represented by a positive semidefinite operator x 
on a Hilbert space of two modes, 

X=j®£('k), (3) 

where I stands for the identity channel and 4' = |4')(4'| 
denotes a density matrix of an infinitely squeezed EPR 
state, 

OO 

|4-) = ^|nn). (4) 

n—0 

The input-output transformation pout = ^iPin) can be 
expressed as 

Pout = Trin [pf^ 0 /%] , (5) 

where T denotes transposition in Fock basis, I stands for 
the identity operator, and Trin denotes partial trace over 
the input mode. In Fock basis, the formula (§ explicitly 
reads, 

OO OO 

Pout ,m,n — EE Xkm,lnPm,k,l ■ (6) 



FIG. 1: (a) Ancilla-assisted quantum process tomography 
[3 El HQ] of a single-mode operation £. The operation is ap¬ 
plied to one part of input two-mode squeezed vacuum state 
1^) and both output modes are measured with balanced ho¬ 
modyne detectors HD^ and HDs- (b) Homodyne measure¬ 
ment of quadrature of mode A of the two-mode squeezed 
vacuum state |^) prepares mode B in a coherently displaced 
and rotated Gaussian squeezed state, (c) The ancilla assisted 
tomography is therefore equivalent to probing the operation 
£ with a suitably chosen ensemble of single-mode squeezed 
states. 


Here pm,n = (w|p|n) and Xkm,in = {km\x\l^- Identity 
channel I is isomorphic to the EPR state ffl, xi = 

and XT,km,In — ^km^ln- 

The input-output transformation ([^ can be also for¬ 
mulated for phase-space representations. Let Win{x,p) 
and Woutix,p) denote the Wigner functions of input and 
output density operators pin and Poutj respectively, and 
let Pin, Xout,Pout) denote the Wigner function of 

operator y. The partial trace ([^ can be rewritten as an 
integral over the phase space of the input mode, 

/ OO POO 

I tVin{Xin, Pin) 

-OO J —OO 

^ (^in: Pin; ^out , Pout )dXindpin. 

where IFin(a;in, — Pin) is a Wigner function of the trans¬ 
posed input state p^. 

Our goal is to establish a procedure for determina¬ 
tion of the matrix elements Xkm,in from experimental 
data. Formula ^ suggests that this could be achieved 
by probing the quantum operation £ with one part of 
the EPR state jT). To make this continuous-variable 
ancilla-assisted quantum process tomography puni ex¬ 
perimentally feasible, the unphysical infinitely squeezed 
EPR state may be replaced with a two-mode squeezed 
vacuum with finite squeezing [S], 

OO 

ie) = \/r^5]A-M, (7) 

n—0 

where A = tanhr, and r denotes the squeezing con¬ 
stant. After some algebra, we find that the elements 
of the quantum process matrix y can be determined as 
properly rescaled elements of the output two-mode state 
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cr^ = I®£{C). where ^ = | 0 (CI, 

= ( 8 ) 

If both modes of the output state would be mea¬ 
sured with homodyne detectors, see Fig. 1(a), then the 
matrix elements could be reconstructed by quan¬ 

tum homodyne tomography Pi. Let rjA and rjB de¬ 
note the overall detection efficiency of balanced homo¬ 
dyne detectors HD^ and HD^, respectively. A detector 
with efficiency rj can be modeled as a lossy channel with 
transmittance rj followed by an ideal detector with unit 
efficiency. The detectors measure rotated quadratures of 
modes A and B, which are specified by angles 9 and 


respectively, 

= VVAi^A cos9 + pa sin 6 ») -b a/1 - VaXa^b^c^ 
x% = coscj)+ PB sincj)) + ^1- 

( 9 ) 

Here xj and pj denote the amplitude and phase quadra¬ 
tures of mode J, [xj,pk] = iSjK, and and Xg ,^^^ 

represent quadratures of auxiliary vacuum modes. This 
measurement samples the joint quadrature distribution 
P{x\,x'^g-,9,riA,(j)T'nB) and the matrix elements Xkm,in 
can be determined by averaging the so-called pattern 
functions over the quadrature statistics [57H55] . 


\km,ln 


47r2^fc-l-/ 


oo 27r oo 27r 



— oo 0 —oo 0 



Here fm,nix,rj) represent the loss-compensating single¬ 
mode pattern functions for density matrix elements in 
Fock basis. Explicit analytical expressions for fm,n{x,r]) 
are provided in Ref. |39j . Since these expressions are 
rather cumbersome, we do not reproduce them here. We 
only note that the pattern functions /m,n(a;, 77 ) are well 
defined for 77 > ^ and they diverge when 77 —|. 


III. SINGLE-MODE PROBE STATES 

In this section, we will propose a procedure for 
quantum process tomography with single-mode squeezed 
probe states. In particular, we will exploit the fact that 
the ancilla assisted process tomography with a two-mode 
squeezed vacuum state |^) and individual single-mode ho¬ 
modyne measurements on the output modes is equivalent 
to preparation of a specific ensemble of displaced and ro¬ 
tated single-mode squeezed states of mode B, followed by 
probing the quantum operation £ with these states |23j . 
To see this equivalence, we rewrite the joint quadrature 
distribution as 

Pix%x‘^g-9,T]A,(t),r]B) = P{x^A'^e,r]A) 

xP{x%;(j),r]B\x^A'^6,r]A), 

( 11 ) 

where P{xA',9,r]A) is the probability density of 
measurement outcomes on mode A, and 

P(x^; ((), 77 b|x^; 0 , 77 a) is the conditional probability 
density of measurement outcomes of quadrature Xg 
on mode B provided that a particular measurement 
outcome x^ was obtained on mode A. 


fk,i{xA^r]A)fm,n{x%,r]B)P^'' OSgd"! ^)<l>dx\d9dx%d(t). 

, _^ 


Since mode A is in a thermal state, the probability 
P{xA'^9,r]A) does not depend on 6, and all quadratures 
Xa exhibit Gaussian distribution with zero mean and 
variance 


^A = ^ [VA cosh(2r) -f 1 - 77a] . (12) 


This formula accounts for imperfect detection with effi¬ 
ciency rjA, and 1 cosh( 2 r) is the variance of quadratures 
of mode A of the pure two-mode squeezed vacuum state 
Q. Explicitly, the probability density reads 


P(x^;6»,77a) 





(13) 


Homodyne detection of quadrature x^on mode A of 
the two-mode squeezed vacuum state Q prepares the 
other mode B in a coherently displaced squeezed state 
with squeezing ellipse rotated by angle —9, see Fig. 1(b). 
This rotation follows from the identity 

UA{9)UBi-9m = 10 , ( 14 ) 

where U{9) = is a unitary phase shift operator. 

Measurement of a rotated quadrature x^ on mode A is 
thus fully equivalent to measurement of quadrature xa, 
followed by rotation of mode B by —9. The covariance 
matrix of the conditionally prepared state does not de¬ 
pend on the measurement outcome x^, and the coherent 
displacement d is linearly proportional to the measure¬ 
ment outcome. 

Let V- and V+ denote the variances of squeezed and 
anti-squeezed quadratures of the conditionally prepared 
state, and let d denote the coherent displacement od the 
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squeezed quadrature of this state. It follows from the 
above discussion that, without loss of generality, we can 
assume 0 = 0 in our derivation of V-, V^, and d. It is 
convenient to collect the quadrature operators of modes 
A and B into a vector z = (xa,Pa,xb,Pb) and define a 
two-mode covariance matrix = {AzjAzk -f AzkAzj), 
where Az^- = Zj — (zj). Covariance matrix of a two-mode 
squeezed vacuum state Q whose mode A was transmit¬ 
ted through a lossy channel with transmittance rjA reads, 


/2 Ua 

0 

K 

0 

0 

2 Fa 

0 

-K 

K 

0 

214b 

0 

V 0 

-K 

0 

2 yB 


(15) 


where Vb = ^ cosh(2r) and K = y^r^sinh(2r). 

Since there are no correlations between the xa and pb 
quadratures, measurement of xa does not influence ps, 
whose variance remains equal to Vb and (ps) = 0, 

14+= - cosh(2r). (16) 


In contrast, the measurement of xa will reduce fluctua¬ 
tions of xs due to the correlations between xa and xb- 
The resulting (conditional) variance V- of xb can be cal¬ 
culated by minimizing the variance of xs — gxA over a 
tunable gain g. The optimal gain reads popt = K/(2Va), 
which yields 


1 PA -f (1 - pa) cosh(2r) 

2 r]A cosh(2r) + 1 — gA 


(17) 


Moreover, the coherent displacement d of the condition¬ 
ally prepared state of mode B is given by (xb) = Popt^A^ 
which explicitly reads 


d = 


^/pjsinh(2r) 

PA cosh(2r) -I- 1 - PA 


(18) 


The variances V- and V+ of squeezed and anti¬ 
squeezed quadratures of the probe single-mode state de¬ 
termine the effective detection efficiency pA of HDa and 
the parameter A of the (virtual) two-mode squeezed vac¬ 
uum state Q. By inverting formulas (16) and ( [T7| ), we 
get 


VA = 


2(V+ - lA) 

(214+ - l)(2y_ -h 1) ^ 


(19) 


and 




( 20 ) 


The effective detection efficiency pA > | if and only if 
the probe state is squeezed and 14_ < 1. This establishes 
single-mode squeezing as a valuable resource for contin¬ 
uous variable quantum process tomography. If the probe 
state is pure, V+ = l/(4t4_), then pA = 1. If t4_ < 1 


then the efficiency is a decreasing function of 14+ and in 
the limit V+ —)■ oo we get pA = 1/(1 + 2V-). The co¬ 
herent displacement (IT^ of mode B can be expressed in 


terms of the quadrature variances as follows. 


d = ^2iV+-V.)J l^: ^ ] xA. 


214, 


1 


( 21 ) 


Formula 0 together with the above results 
suggests that the joint quadrature distribution 
P(x^, x^; 0, Pa, (/, Pb) can be sampled as follows. 
Generate random x^ drawn from the Gaussian dis¬ 


tribution (13) and a random 6 and (j) drawn from a 


uniform distribution in the [0, 27r] interval. Prepare a 
single-mode squeezed Gaussian state with variances 14_ 
and 14+ and displacement d, rotated in phase space by 
—9, as illustrated in Fig. 1(b). Send this probe state 
through the quantum channel £ and measure a rotated 
quadrature Xg of the output state with a homodyne 
detector. Note that the squeezing properties of the 
input probe states do not depend on x^ and 6, hence a 
source producing squeezed states with a fixed amount of 
squeezing and anti-squeezing is sufhcient. 

The improvement achieved by squeezed probe states in 
comparison to coherent probe states comes at a cost of 
somewhat increased experimental difficulty. In particu¬ 
lar, the variances 14+ and t4_ of the probe squeezed states 
need to be precisely characterized, which can be achieved 
by routine homodyne detection, and these parameters 
have to be kept constant during the whole tomographic 
measurement. Moreover, the orientation of the squeez¬ 
ing ellipse should be fully under control and tunable to 
any required angle 9. Finally, the ability to coherently 
displace the squeezed state is also required, which can 
be achieved e.g. by mixing it with an auxiliary coherent 
beam on a highly unbalanced beam splitter |40) . 


IV. TOMOGRAPHY OF QUANTUM 
MEASUREMENTS 

The proposed method can be also adapted to tomo¬ 
graphic characterization of quantum measurements [Sl- 
I331I1IH111- Gonsider a detector D which can respond with 
K different outcomes. Each outcome is associated with 
a POVM element 11^ and the probability to observe an 
outcome k for input state p reads p{k) = Tr[n^p]. Con¬ 
sider now an ancilla-assisted quantum detector tomogra¬ 
phy [43l [44] , where the detector is probed with one part 
of a two-mode squeezed vacuum state, see Fig. 2(a). The 
conditionally prepared state of mode A corresponding to 
measurement outcome k on mode B can be expressed as 

OO OO 

p'==(l_A2)^^A™+"<,„|n)(m|. (22) 

m—0 n—0 

The state is not normalized and its trace is equal to prob¬ 
ability of observing the outcome fc, 

Tr(p'=) = (e|/A®n^lO- 


(23) 
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e) 





B 



(a) 


(b) 

Od 


The matrix elements of 11^ can then be immediately ob¬ 
tained from Eq. (24), where the parameter A = tanh r is 
determined by the variance of anti-squeezed quadrature 
V+ of the probe state, see Eq. (16). 


FIG. 2: (a) Ancilla-assisted quantum detector tomography 
|43[ I44| . An unknown detector D is probed with one part 
of two-mode squeezed vacuum state |^). Information about 
POVM element 11*^ associated with measurement outcome k 
of detector D is imprinted into the corresponding conditional 
state p*’ of mode A, which can be characterized by homodyne 
tomography, (b) An equivalent scheme where the detector is 
probed with an ensemble of Gaussian squeezed states. 


Eormula (221 implies that the information about the 
POVM element 11^ is imprinted into the conditional state 
. In particular, we have 


Ij/c x-(m-l-n) A _ x2\ ^ fc 

^ ^ J rn,mf 


(24) 


in analogy with Eq. (j^ . The conditional states can be 
characterized by homodyne detection on mode A, which 
would provide sufficient data to reconstruct the density 
matrix elements „. 

Similarly as for quantum operations, probing with one 
part of two-mode squeezed vacuum can be replaced by 
probing with single-mode squeezed states, see Eig. 2(b). 
Let p{k\x^^,9) denote the probability of outcome k for a 
probe state with displacement and rotation s peci fie d by 
parameters a;^ and 6, respectively, c.f. Eqs. (16), (17), 
and (18). The statistics of homodyne measurements on 
p^ is governed by P{x^,O,r]A)pik\x^,0), where rjA is a 
function of the variances of squeezed and anti-squeezed 
quadratures of the probe squeezed state, see Eq. (IT^. 


The density matrix elements of p^ can be obtained by 
averaging appropriate pattern functions over the quadra¬ 
ture statistics. 


V. CONCLUSIONS 


In summary, we have proposed a procedure for to¬ 
mographic characterization of continuous variable quan¬ 
tum operations which employs homodyne detection and 
single-mode squeezed probe states with a fixed degree of 
squeezing and anti-squeezing and a variable displacement 
and orientation of squeezing ellipse. We have shown that 
the elements of quantum process matrix x in Fock basis 
can be estimated by averaging suitable pattern functions 
over the homodyne data. The pattern functions are well 
behaved provided that the probe state is squeezed and 
V- < ^. For the sake of simplicity, we have considered 
tomography of a single-mode operation £. However, the 
method can be straightforwardly extended to multimode 
operations. For tomography of V-mode operation, one 
would have to use N independent single-mode squeezed 
states and measure each output mode with an indepen¬ 
dent homodyne detector. While we have focused on lin¬ 
ear reconstruction procedure based on pattern function 
formalism, other methods of data processing would be 
also possible. For instance, one may utilize the widely 
employed maximum-likelihood estimation, or other ap¬ 
proaches. Given its relative simplicity and practical feasi¬ 
bility, the present procedure is likely to find applications 
in the characterization of continuous variable quantum 
operations and measurements. 
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